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I present here a concise review of the experimental results obtained at the Relativistic 
Heavy Ion Collider (RHIC), which shed light on the hot and dense quark gluon matter 
produced at these high temperature and density conditions. 
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1. Foreword 

This review has little to do with the topic ol the symposium it was given in, namely 
super-symmetry at the Large Hadron Collider (LHC). Indeed, it deals with standard 
QCD at high density and temperature and how it was widely probed at the Rela- 
tivistic Heavy Ion Collider (RHIC). In fact, an intimate theoretical link does exist, 
as first stated by MaldacenaH between strongly coupled four dimensions super- 
symmetric0 Yang Mills theories and weakly coupled type lib string theories on 
Anti-deSitter five dimensions space. This so-called AdS/CFT correspondence allows 
one to compute some properties (viscosity/entropy ratio, quenching parameters...) 
of the quark gluon matter, as if it were a black hole. Even if such a connection 
exists on the theory side, I will not comment further on this matter and rather 
review various experimental facts that highlight certain properties of the matter 
produced in relativistic (up to ^Jsnn — 200 GeV per pair of nucleons) heavy ion 
collisions. The early and most striking RHIC results were summarized in 2005 by 
the BRAHMS, PHENIX, PHOBOS and STAR experiments in their so-called white 
that will be largely referenced thereafter. 

2. Multiplicities and densities 

The first obvious thing that comes out of a heavy ion collisions is a lot of parti- 
cles. The number of charged particles was measured for various collision energies 



a But is QCD super-symmetric? 
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and centralities by the four RHIC experiments, and in particular by the dedicated 
PHOBOS collaboration over a broad range of 10.8 units of pseudorapiditjEl (as 
illustrated by Fig. [T] on the left). At midrapidity, the number of charged particle 
reaches dNch/drjl^—o ~ 670 in the most violent Au+Au collisions and they sum up 
to about 6000 particles (of any charge) over the full rapidity range. As illustrated 
by the right part of Fig. [1] these huge numbers were in fact lower than expected 
from various simple models, extrapolating lower energy results (for more details 
and complete references, see section 2.1 in Ref. 4 or Ref. [7J). This moderation of the 
produced particles is an indication that the gluon density in the initial state starts 
to saturate, or similarly to be shadowed. In other words, low momentum gluons 
from neighbor nucleons overlap and recombine. Indeed, a strong gluon shadowing 
in the first model on the figure (HIJING) is necessary to reduce the multiplicities 
to the observed level. Such a saturated initial state is also described in the Color 
Glass Condensate (CGC) framework, which is fairly able to reproduce the particle 
multiplicities (the McLerran-Venugopalan "McLV" last model is such a saturation 
modepJ) . 

Another glimpse into the CGC is provided by the high transverse momentum 
(pr) particle suppression observed in d+Au collision at the highest (up to r] = 3.2) 
pseudorapidity by the BRAHMS experimentP. Indeed, higher rapidity allows one 
to probe lower gluon momentum fractions x in the gold nucleus and one clearly 
sees an increasing suppression with rapidity (see section 7 in Ref. [5] and references 
therein). 

These results show that the (initial) matter is gluon saturated. 
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Fig. 1. Left: Charged particle multiplicities as a function of pseudorapidity for various energy and 
centrality Au+Au collisions from the PHOBOS experiment. Right: The top most dN ch /dy\ y= o 
compared to various models. 
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The second obvious manifestation of the collision violence is the transverse (i.e. 
unboosted by the initial parton longitudinal momenta) energy liberated. Measur- 
ing it allows one to estimate the en ergy density e of the medium after a given 
time To, through the Bjorken formula^ £ = dE T j f dy\ y= o/ 'tqAt , where At is the 
transverse area of the collision. The four RHIC experiments measure consistent val- 
ues of dET/dy\ y=0 that correspond to an energy density of at least 5 GeV/fm 3 at 
To = 1 fm/c. The question of the time to be considered is not trivial, but 1 fm/c is 
a maximum if one cares about the earliest as possible thermalized medium. Indeed, 
hydrodynamics models provide thermalization times between 0.6 and 1 fm/c, while 
the formation time is estimated to be 0.35 fm/c and the nucleus- nucleus crossing 
time is 0.13 fm/c. For a detailed discussion of energy density and time scale esti- 
mates, see section 2 of Ref. |3] What matters here is that the lower energy density 
estimate is much higher than the threshold for the transition to a quark gluon 
plasma, as predicted by QCD on the lattice^} e c ~ 1 GeV/fm 3 . 

This tells us that the matter should be deconfined, i.e. made of free quarks 
and gluons. The following sections review the main measurements that indicate that 
it is indeed the case. 



3. High Transverse Momentum Suppression 

Fig. [2]is an illustration of the first a nd most striking QGP signature seen at RHIC, 
namely the quenching of 

jetJ™. 

Displayed is the nuclear modification factor 
Raa defined as the yield of particles seen in A+A collisions, normalized by the 
same yield from p+p collisions scaled by the average number of binary collisions 
corresponding to the considered centrality: Raa = cINaa/ (N C oli) x dN pp . Hard 
processes (high px particles in particular) are expected to respect such a scaling 
(Raa = !)■ This is indeed the case of the direct photorP^l (purple squares) up to 
13 GeV/c, while the corresponding ir° (orange triangles) and r/ up to 10 GeV/c (red 
circles) are suppressed by a factor of fiv^B. This is understood as an energy loss of 
the scattered partons going through a very dense matter, and producing softened 
jets and leading (high px) particles. The medium is so dense that it cannot be made 
of individual hadrons, but rather of quarks and gluons. Gluon densities of the order 
of dN 9 /dy ~ 1100 are needed to produce such a strong quenching^! 

High px suppressio n is see n for various particles with various pt reaches and 
by the four experiment JHSEE1. R gets stronger for more central collisions. It is not 
observed in d+Au collisions (in particular for neutral pions^ to be compared to 
the ones on Fig. [5]) where a moderate enhancement is even seen as a function of px, 
probably due to multiple scattering of the incoming partons providing additional 
transverse momentum (the so-called Cronin effect). 



b It is to be not ed that PHENIX has released preliminary data up to 20 GeV/c for ir° and 18 GeV/c 
for photons!^. While the 7r° remain at Raa ~ 0. 2. ph otons start to deviate below unity, possibly 
because of the nucleus to proton isospin differenc a^l which has nothing to do with QGP. 
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PHENIX Au+Au (central collisions): 




2 4 6 8 10 12 14 16 

p T (GeV/c) 

Fig. 2. Nuclear modification factors for photon, r) and n° for central collisions, from the PHENIX 
experiment . 

This quenching of high pr particles shows that the matter they traverse is 
dense. 

4. Back to Back Jets 

Another way to look at jets is to consider back to back high transverse momentum 
hadron correlations. Fig. [3] shows the measurements of such correlations for various 
collision types performed by the STAR experiment and reported in section 4.2 of 
Ref. [5j Displayed are the azimuthal distributions of hadrons around a "trigger" 
particle of high enough pt to reflect the main direction of jets (4 GeV/c for the 
trigger particle and 2 GeV/c for the others in this example). In p+p collisions (black 
histogram), one clearly sees particles belonging to both the narrower same (A</> = 0) 
and broader opposite (A(f> = tt) jets, while in central Au+Au collisions (blue stars) 
the away-side jet disappears^. This is also attributed to jet quenching, the away- 
side jet being absorbed by the dense matter produced at RHIC. As for the high 
Pt suppression we saw in the previous section, this effect is not observed in d+Au 
collisions (red circles) where away-side hadrons are clearly distinguishable^^. 

Jet-induced hadron production has been further and extensively investigated at 
RHIC and various effects corroborate the jet quenching hypothesis, among which: 

• In Au+Au collisions, the away-side disappearance grows with centrality. In 
fact, the most peripheral collisions exhibit a very similar away-side pattern 
as in p+p and d+Au collisions. 
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The jets emitted in the reaction plane are less suppressed than in the per- 

lOI I 

pendicular direction, where they have more matter to traversed In lact, 
the high px (near-side) particles we see in central Au+Au collisions are 
likely to come from the periphery, the "corona" , of the collision. 
By lowering the pt r equirements (down to ~1 GeV/c), one can find back 

IOOI 

the away-side jets^ 1 . 

These weakened away-side jets are depleted at Acf> = ir and exhibit two dis- 
placed maxima around Acf> = tt± 1.1 radians^! This camel-back or conical- 
like shape provides insight in the quenched parton interac tion s with the 
medium. Various scenarios are proposed, such as radiative loss^S Cerenkov- 
like or Mach-cone emissiona^. The later allows one to compute an average 
speed of sound in the medium of eg ~ 0.45. 

Prelim inar y analyses of three particles correlations also exhibit the conical 
patterrM 

The near-side jet exhibits a "ridge" along pseudorapidity (thus perpendic- 
ular to the azimuthal stru cture) that suggests the jets are indeed flowing 
with the expanding mattep2212ZI 



In brief, these high px dihadron correlation studies show that the matter is 
opaque to jets to a first approximation, and clearly modifying their remaining 
structure. 
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Fig. 3. Dihadron azimuthal correlations in p+p, d+Au and Au+Au central collisions, from the 
STAR, experiment. 
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5. Elliptic Flow and Ideal Hydrodynamics 

Speaking of azimuthal correlation, it is noticeable that for moderate centralities, 
overlapping colliding nuclei form an almond-shape area. It is then relevant to look 
at the "elliptic flow" of particles, namely the second Fourier harmonic V2 of the 
azimuthal distribution: dN/dcj) = A* r (l + 2'f; 1 cos(0) + 2v 2 cos(20) + . . .). Experimen- 
tally, V2 happens to be positive, meaning that the particle emission is enhanced in 
the plane of the reaction (along the smaller axis of the almond) with respect to 
the out-of-plane emission (along the larger axis). This reflects pressure gradients, 
i.e. strong interactions, that must exist at the very early stage of the collision to 
provide more transverse momentum to the emitted hadrons along the shortest axis. 
Moreover, the rather large values (up to v-i ~ 20 % at 2 GeV/c) of the elliptic flow 
measured at RHIC contradict hadronic transport models (for instance accounting 
for only ~ 60 % of the observed value^j). On the contrary, ideal hydrodynamical 
models (for a list see section 3.5 of Ref. |3]) that assume a QGP equation of state, 
a high energy density (e ~ 20 GeV/fm 3 ) and fast equilibration time (r ~ 0.6 to 
1 fm/c) fits reasonably well a broad selection of data: 



The transverse momentum dependence of elliptic fl ow is reproduced up 
to 2 GeV/c, and properly ordered for various specieJ^^ (from pions to 
cascadesjj. 

These V2(pt) scales with the eccentricity ((y 2 ) — (x 2 )/(y 2 ) + (x 2 )) of the 
reaction for various collision systems, centralities and energies, underlining 
the facts that elliptic flow does reflect the very early stage of the reaction 
and that thcrmalization must arise ra pidljEH 

Hydrodynamics pressure gradients also imply a scaling by the transverse ki- 
netic energy. While this property is verified for lowpT (less than ~ 1 GeV/c) 
hadrons, it extends its validity to much higher p-r when one divides both v-i 

jo -l I 

and pt by the number n of constituent quarka^M This r esult holds for pi- 
ons, kaons, protons, A, H, but also for the 4> mesoniPZEH] ( f baryonic-level 
mass) and deuterona^ (with n = 6), as shown on Fig. [4] 
The adjunction of even a low viscosity in hydrodynamical models deterio- 
rates their fits to the data, in particular by moderating w 2 as pr growiP^I 
(departing from ideal hydrodynamics around pt ~ 1 GeV/c). The matter 
created at RHIC must then have a very low viscosity and was thus qualified 
as a "perfect liquid". 

The transverse mass spectra, i.e. the radial flow, are also reproduced by hy- 
drodynamical models (with kinematic freeze-out temperature of ~ 100 MeV 
and transverse speed of (/3t) ~ 0.6 for the most central Au+Au collisions^). 



c Being faster, higher pt particles share less the collective behavior of the bulk, which does not 
mean they do not see it, since we saw in the previous sections that they are very suppressed by 
this dense matter they traverse. 
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Fig. 4. Scaling of the elliptic flow parameter V2 versus transverse momentum (left) or kinetic 
energy (right) for various particles. Both quantities are divided by the number of consistent quarks. 

This high degree of collective ideal hydrodynamical behaviors, setting up at very 
early times and exhibiting a low viscosity, tells us that the matter is strongly 
interacting, in a liquid-like manner. 

6. Baryons and Mesons 

We saw on Fig. [4] that dividing by the number of constituent quarks helped the 
elliptic flow parameter vi to scale with transverse kinetic energy at moderate pt 
(2 to 4 GeV/c). This is not the first observable to exhibit a parton-like scaling. 
Indeed, the nuclear modification factor Raa also shows a different pattern between 
baryons and mesons (including the <j> which is of baryonic- level mass) , in the same 
Pt range. Fig. [5] from Ref. [5] shows the central to peripheral ratios (scaled by the 
number of collisions) of mesons (left) and baryons (right). It is first noticeable that 
the relevant property to determine the fate of these intermediate pt particles is 
their baryonic/mesonic nature. Moreover, the fact that the baryon peak production 
is higher and lays at ~ 3/2 times the mesonic one, suggests that a quark coalescence 
or recombination mechanism is at play. 

To test this hypothesis, the p/n + and p/tt~ ratios can be studied in detail. 
Baseline p+p and peripheral Au+Au collisions exhibit very similar patterns, while 
the p/ir ratio is clearly enhanced in the moderate pt range. Th ese particle ratios 
are equally reproduced by coalescence or recombination approach^!. 

Added to the elliptic flow versus transverse kinetic energy scaling (Fig. [4]) , and 
to the partonic strength of jet quenching (Fig. [2|), this result suggests that the 
matter is of partonic nature. 
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Fig. 5. Central to peripheral ratios of various mesons (left) and baryons (right) as a function of 
transverse momentum, as measured by the STAR experiment. 



7. Heavy Flavors Flow and Quenching 

Being heavier, charm or bottom quarks are produced earlier than the light flavors, 
and their production yields can be in principle calculated by perturbative QCD. 
They are thus considered as good probes of the plasma earliest times. As we saw for 
light quarks, the nuclear modification factors and elliptic flow are good observables 
of the medium effects on produced particles. Fig [5] shows both quantities for elec- 
trons from heavy flavor decays (blue c ircle s) and ir° (shaded band or red squares), as 
measured by the PHENIX experiment. It is to be noted that even if the STARp3 
and PHENIX experiments disagree on the charm cross-section, they do agree on V2 
and Raa- They both see that high pt heavy quark^ are quenched by a factor of 5 
and that they do exhibit a significant flow (up to 10 %, while pions reach 20 %). 
As for light flavors, both observables reveal a strong coupling to the medium. 

These were surprises. Energy loss in a gluon medium was expected to be reduced 
for heavy quarks. Indeed, in order to reproduce the data, one would need a much 
higher gluon density than the one required for light flavors (dN g /dy ~ 3500 versus 
1100, neglecting less quenched beauty decayJ^Hl). Various hypotheses are made to 
reinforce the heavy quark quenching (adjunction of elastic energy loss, change in 
the charm/beauty ratio, modification of the strong coupling constant. . . ). Another 
approach is to quantify the medium effects by transport or diffusion coefficients. 
The models displayed on Fig. [6] follow such approaches. With rather high values of 
these coefficients (q = 14 GeV 2 /fm in mode l I) they roughly manage to reproduce 
the amount of suppression and flo#2l4Q|41j 



d The low pj> (up to 1.5 GeV/c) dominant yield scale with the number of collisions (Raa — 1) as 
expected. 



February 1, 2008 13:33 WSPC/INSTRUCTION FILE 2006.03. 13. Cairo 



What's the matter at RHIC? 9 



Interestingly, and in order to illustrate the black hole/QGP connection men- 
tioned as a foreword, the diffusion coefficient used in the models labeled II and III, 
corresponds to viscosity over entropy ratios that fall close to the lower quantum 
bound of h/AnkB as derived through the AdS/CFT correspondence^!. 
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Fig. 6. Heavy flavor decay electrons in-medium behavior as measured by the PHENIX experi- 
ment, compared to it and models. Top: Quenching in the most central Au+Au collisions. Bottom: 
Minimum bias elliptic flow. Both are as a function of transverse momentum. 



It is fair to say that the way the RHIC matter impacts heavy quarks is not 
perfectly understood yet, but it is also clear that it is strong. To that extend, I will 
dare to say that the matter is "tough", tough to understand as well as tough 
enough to shake the heavy flavors. 
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8. Quarkonia Suppression 

We just saw that the bulk (low pr) charm production scales to first order with the 
number of binary collisions. This forms a good baseline for the study of bound states 
made of charm-anticharm quarks, the more stable of which being the J/ip particle. 
In fact, charmonia were predicted to melt in the QGP, due to Debye screening of 
the color charge^!. Furthermore, J/ip suppression was indeed observed at lower 
energy (•>/ s nn = 17.3 GeV) by the NA50 experiment!^ and is the main signature 
that led CERN to claim for the discovery of QGP. It was thus very awaited at 
RHIC energies. Fig. |7[ shows J/ip nuclear modification factors as measured by the 
PHENIX experiments^], for both midrapidity (red circles, \y\ < 0.35) and forward 
rapiditiy (blue squares, 1.2 < \y\ < 2.2), compared to the NA50 result (black 
crosses) and as a function of centrality (given here by the number of participants 
Npart)- The midrapidity result is surprisingly similar to the pattern observed by 
the NA50 experiment, which also lies close to midrapidity (0 < y < 1). There is 
no fundamental reason for this to happen. First, the energy density for a given 
Npart is much higher at RHIC and should further melt quarkonia. Second, J/ip 
are known to be su ppressed by regular nuclear matter as it is seen in p+A or 
d+A collisionsPESl and this normal suppression should be different. In order to 
compare the two energy regim es, one first need to subtract the regular nuclear 
matter effects. Such an attemplF^ is shown on the right part of Fig. [7] At RHIC, 
these effects are poorly constrained by a relatively low statistics d+Au dataseP^Q 
which significantly increases the uncertainty on the J/ip survival probability in the 
produced matter. 

Anyway, we clearly see that J/ip are suppressed beyond normal nuclear effects, 
both at SPS and RHIC (especially at forward rapidity). Then, once these effects 
are subtracted, these facts remain: 

• SPS and midrapidity RHIC J/ip suppression are possibly still compatible. 
This lead to the hypothesis that direct J/ip do not melt neither at SPS nor 
at RHIC, but that onl y th eir feed-down less-bound contributors (xc and ip') 
disappear in the QGP^Sl. 

• J/ip are more suppressed at forward rapidity. This seems to contradict 
all models based on density-induced suppression, in particular the original 
Debye screening hypothesis, as well as the sequential melting scenario of 
excited states suggested above. 

Two ideas exist to explain this last surprising feature. First, gluon saturation 
could further suppress J/ip at forward rapidity, by playing a larger role than the one 
simply extrapolated from d+Au collisions (such an effect was computed for open 
charm in the CGC framework^^) . Second, J/ip could be recreated in the plasma 

e As well as by the fact that different p+p references were used for d+Au and Au+Au nuclear 
modification factors, which accounts for the ~ 30 % global uncertainty quoted on the figure. 
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Fig. 7. J/tp suppression measured by the PHENIX and NA50 experiments, as a function of 
centrality, given by the number of participants. Left: nuclear modification factor. Right: J/ip 
survival probabilities after normal nuclear effects subtraction. 



by recombination of independent ch arm and antich arm quarks (a large variety of 
recombination or coaelescence 

model j50 | 51 | ft2 | 53 | 5 41 exists) 

. These two ideas do not 

provide quantitative predictions of the nuclear modification factors (recombination 
models suffering from the lack of input charm quark distributions). Other observ- 
ables (pt and rapidity dependencies) are also availablJ^ but so far, they do not 
allow to rule out any possibility. New measurements are thus needed to further 
understand J/tp suppression at RHIC (higher statistics d+Au, J/tp elliptic flow, 
feed-down contributions...). 

However, we do not need them to reckon that J/tp do melt beyond normal 
nuclear effects. This is a sign that the matter is deconfining. 



9. Thermal Radiation 

Last but not least, a thermalized matter as the one suggested by the strong elliptic 
flow should emit its own thermal radiation. We saw on Fig. [2] that photons are 
unmodified by the medium and the nuclear modification factor is compatible with 
unity. This holds for px > 2 GeV/c, but lower px photons exhibit an enhance- 
ment when compared to perturbative, next-to-leading order QCD predictions^! 
This is further illustrated for the most central Au+Au collisions on Fig. [8] on which 
the lowest px photon^ deviate from the prompt contribution from a NLO pQCD 
calculation (dashed line) . They are con siste nt with the addition of a thermal contri- 
bution. Various hydrodynamical modelJ^ fairly reproduce the data assuming early 
(typically at a time of the order of 0.2 fm/c) temperature of 400 to 600 MeV, well 

'Note that the internal conversion 7* yield is a PHENIX preliminary resull!551. 
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above the critical temperature of T c = 190 MeV provided by lattice QCE^ as the 
phase transition boundary to a quark-gluon plasma. 




j r\-7 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

10 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 

P T (GeV/c) 

Fig. 8. Thermal + perturbative QCD fits to the photon yield in central collision, as seen by the 
PHENIX experiment (the internal conversion 7* yield is preliminary). 



This result suffers the lack of experimental p+p reference, but if NLO pQCD is 
taken as a baseline, we do see thermal photons that demonstrate that the matter 
is hot. 



10. Conclusions 

Even if we haven't (yet) observed any sharp change in the behavior of the Au+Au 
observables related to the predicted phase transition, nor numbered degrees of free- 
dom, it is clear that the matter produced at RHIC behaves very differently than 
ordinary hadronic matter. Indeed, to answer the question raised by our title, we 
saw that the matter is gluon saturated, dense and opaque, strongly interacting and 
liquid-like, partonic and deconfining, "tough" and hot. It is thus very likely to be 
formed by deconfined quarks and gluons. 
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